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Abstract 

We calculate the widths of the various pentaquarks in a relativistic model in which the pentaquark 
T^j- ' is considered to be composed of a scalar diquark and a spin 1/2 triquark. We consider both positive 

o ' 

■ and negative parity for the pentaquark. There is a single parameter in our model which we vary 

o 



a: 



and which describes the size of the pentaquark. We obtain quite small widths for the decay 
+ — > N + K + and for — > P + D*~ consistent with the experimental situation. For the sum 



a.: 

of the decay widths for E — > 5 + ir + and H — > S + K we find values of the order of 
4 — 8 MeV for pentaquarks of the characteristic size considered in this work. (The experimental 
situation with respect to te observation of the S is somewhat uncertain at this time.) We also 
provide results for the decays N + — > N + it and — ► A + K + . Our model of confinement plays 
an important role in our analysis and makes it possible to use Feynman diagrams to describe the 
decay of the pentaquark. 
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I. INTRODUCTION 



There has been a great deal of interest in the study of pentaquarks and a large number 
of experiments have been carried out [1-11]. The 6 + (1540) which decays to a kaon and a 
nucleon has been seen in several experiments. It has been interpreted as a pentaquark with 
a ududs structure [12]. A pentaquark 6° with the assumed structure ududc has also been 
observed recently. In the present work we will present calculations of the widths of several 
pentaquarks in a relativistic diquark-triquark model which includes a model of confinement 
that we have introduced previously in our study of meson and nucleon structure. 

A number of reviews have appeared. In Refs. [13,14] the experimental evidence for the pen- 
taquark is reviewed. Theoretical and experimental developments are reviewed in Refs. [15- 
18]. A number of theoretical papers have also appeared [12,19-24]. We are particularly 
interested in the work of Ref. [24] in which a diquark-triquark model is introduced to de- 
scribe the pentaquark. We will make use of a variant of that model in the present work, since 
that model lends itself to the analysis of the pentaquark decay using Feynman diagrams. 
The diagram we consider is presented in Fig.l. There the pentaquark is represented by the 
heavy line with momentum P. The pentaquark is composed of a diquark of momentum 
—k + Pn and a triquark of momentum P + k — P^. The triquark emits a quark (u, d or 
s) which combines with the diquark to form a baryon of momentum P^. The final-state 
meson of momentum P — P^ is emitted along with the quark at the triquark vertex. (In 
the simplest model the triquark could be considered as composed of the final-state meson 
coupled to the exchanged quark, which ultimately forms part of the final-state baryon. [See 
Fig. 1.]) 

We have studied the structure of the nucleon in a quark-diquark model in Ref. [25]. In 
that work we considered both scalar and axialvector diquarks, however, in this work we will 
limit our considerations to a nucleon composed of a quark and a scalar diquark. The nucleon 
vertex is described in Ref. [25] for the case in which the diquark is placed on mass shell. 
(We remark that in Fig.l only the quark of momentum k and the triquark of momentum 
P + k — P N will be off mass shell in our analysis.) The mass of the scalar diquark was 
taken to be 400 MeV in Ref. [25]. We use that value here when the diquark does not have a 
strange quark present. It is clear that we need to introduce a model of confinement for the 
pentaquark to carry out our program. Our (covariant) confinement model which we have 
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P+k-P N 



FIG. 1: In this figure the heavy line denotes the pentaquark and the line of momentum — k + P/v 
denotes an on-mass-shell diquark. The line of momentum k represents the quark and P + k — Pjv 
is the momentum of the triquark. In the final state we have a baryon of momentum Pn and a 
meson of momentum P — Pa- 
used extensively in other works [26-28] will be discussed in the next section. 

The organization of our work is as follows. In Section II we will describe our model of 
confinement and in Section III we will discuss the widths of various pentaquarks of positive 
and negative parity. In Section IV we will provide the results of our numerical calculations 
and Section V contains some further discussion. The Appendix contains a discussion of the 
normalization of the various wave functions we have introduced. 

II. A MODEL OF CONFINEMENT 

For the moment, let us consider the case of the + pentaquark with mass 1540 MeV. Our 
pentaquark is composed of a scalar diquark of mass 400 MeV and a triquark of mass 800 
MeV. Since the mass of the pentaquark under consideration is 1540 MeV, the pentaquark 
would decay into its constituents in the absence of a model of confinement. Similarly, the 
nucleon of mass 939 MeV could decay into the scalar diquark of mass 400 MeV and a quark 
whose mass we take to be 350 MeV in this work. 

In earlier work we have introduced a confining interaction which served to prevent the 
decay of mesons or nucleons into their constituents. Our covariant confinement model is 
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FIG. 2: Values of Tg (P°, \ k |) denned in Ref. [29] are shown. Starting with the uppermost curve 
and moving downward, the values of P° are 0, 0.10, 0.20, 0.30, 0.40, 0.50, 0.55 and 0.60 GeV. For 
the last two of these curves r^"(P°, k on ) = 0. Here k 2 on = (Pq/2) 2 - m 2 , with m q = 0.260 GeV. 

described in a series of our papers [26-28]. In that model we solve a linear equation for a 
confining vertex function, T. This function has the following property. Consider the decay 
A — > B + C, in which the hadrons A and C are on mass shell. If we include the confining 
vertex we find the amplitude has a zero when particle B goes on mass shell, so that the 
amplitude for A to decay into two on-mass-shell particles (B and C) is zero. That feature 
may be seen in Fig. 2 which is taken from Ref. [29]. 

For the decay A — > B + C we may introduce a wave function that may be expressed in 
terms of the momentum of the off-shell particle B. If B is a scalar, we have 



¥*(*) = sr ^r fl (A0 . (2.1) 



Note that the ratio of r^(fc) to (k 2 — m 2 B ) is an ordinary function which may often be well 
represented by a Gaussian function. (It is not necessary to include an it in the denominator 
ofEq. (P .) 

In the case of the nucleon we may consider the decay into a quark and a diquark. In 
Ref. [25] we considered both scalar and axialvector diquarks, but for simplicity we will limit 
ourselves to only the scalar diquark. The relevant wave function in this case was given in 
Eq. (3.3) of Ref. [25]: 

* 5 (P, k, s, t) = % } (P, k) , 2mqAi+){k L=u N (P, s) Xt (2.2) 

2E(k) (E{k) + m, 



which we will simplify for the present work to read 

* N (P N , k, s, t) = V N (P N , k)M + \k)u N (P N , s) Xt ■ (2.3) 

The function $?n(P, k) is represented in Fig. 5 of Ref. [25] by a dashed line. That function 
is well approximated by a Gaussian function which we will record at a latter point in our 
discussion. The factor of A^ + \k) = {^ on + m q )/2m q arises from an approximation made in 
Ref. [25]. (Here k on = (E(k), k) with E(k) = [k 2 + m 2 q ] 1/2 .) In Ref. [25] the quark propagator 
was written as 



iS(k) - "'" 



E q (k) 



AW(jfe) AH(-jfc) 



(2.4) 



k°-E q (k) k° + E q (k) 

The second term was neglected in our formalism when we studied the nucleon. (Thus we 
limited our analysis to positive-energy quark spinors.) Since we wish to make use of the 
nucleon wave function determined in Ref. [25], we will continue to include the projection 
\( + \k) in our formalism. Here k is the quark momentum. 

As stated earlier, the diquark of momentum —k + P/v in Fig. ^will be placed on mass 
shell, so that k° = E N (P N ) - E D (P N - k), where E D (P N -k) = \(P N - k) 2 + m 2 D } 1 ' 2 and 
En(Pn) = [Pn + ^at] • That approximation is achieved by writing 
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(Pn — k) 2 — m 2 D + ie 
as described in detail in Ref. [30]. Note that 



-2iriS^[{P N -k) 2 -m 2 D ], (2.5) 



5^[(P N - k) 2 - m 2 D ] = — i =- 5[P N - k° - E D (P N - k)} . (2.6) 

2E D (P N - k) 

The on-mass-shell specification used here arises when performing an integral in the complex 
k° plane [30]. 



III. CALCULATION OF THE WIDTHS OF NEGATIVE AND POSITIVE PAR- 
ITY PENTAQUARKS 

We consider the diagram shown in Fig. ^ Recall that the heavy line of momentum P 
denotes the pentaquark. The line carrying momentum —k + P/v is the on-mass-shell scalar 
diquark and the line with momentum P + k — P/v is the triquark. The momentum k is that 
of an up, down or strange quark, P — P/v is the momentum of the final-state baryon. The 
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pentaquark, final-state nucleoli and final-state meson are all on mass shell. In our analysis 
the diquark is also on mass shell, so that only the triquark and quark propagate off mass 
shell, as noted earlier. (As stated in the last section, the on-shell characterization of the 
diquark arises when we complete the k° integral in the complex k° plane.) 
We now make use of the formula [31] 

dT = \Tl\ 2 ^ ^ (2vr)V(P -k 1 + k 2 ), (3.1) 

(2tt) 3 E N (h) (2vr) 3 2E K (k 2 ) 

where k\ and k 2 are the momenta of the outgoing particles. We may put k\ = Pn and 
k 2 = P — Pn = —Pn f° r a pentaquark at rest. (It is convenient to take P N along the z axis 
when calculating the width.) 

In writing our expression for Y we will represent the product of the quark propagator and 
nucleon vertex function by the nucleon wave function of Eq. (2.3). In a similar fashion we 
will represent the product of the triquark propagator and the pentaquark-triquark vertex 
function by a triquark wave function. (In the vertex, the diquark is on mass shell.) We then 
have to specify the vertex function of the triquark which describes the decay into the quark 
of momentum k and the final-state meson. That scalar part of the vertex is usefully written 
as 

T T (k) = ^ T (k)(k 2 -m 2 q ). (3.2) 

Thus, the wave functions ^jv(^), ^r(k) and ^q(Pn — k) will appear in our expression for 
r when we consider the decay of the 6 + . Note that for the decay + — > N + K + , we have 

y:> _ i'm 2 e -m 2 N + m K 



n = i 2 ;: K ' -< < 33 » 



which yields P% = 0.0722 GeV 2 , or |P^| = 0.269 GeV. 
We find that the width is given by 



1 f dk 1 f dk' 



e 



2 J (2tt) 3 2E D (k-P N )J (2tt) 3 2 E D (k> - P N ) N e N T N N (2m N )(2m e ) 
x * N (k) V N (k') V e (P N - k) V e (P N - k') {k 2 - m 2 q )V T (k) (k 12 - m\)^ T (k'\ 

x Tr [($ - P N + P + m T ) (^ on + m q ) (P N + m N ) 7 ° 

(Kn + m q ) {fi-P N + P + m T ) 7 ° (P + m e )]p, 



(3.4) 
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where p is the phase space factor. We find for the decay of the + , 

P = 7i \ P N\ 

2tt me 
= 0.0261 GeV, 

since \P N \ ~ 0.269 GeV. 

The factor f N + m N and P + m® in the trace arise from the relations 

?n + m N 



(3.5) 
(3.6) 



2m 



N 



= ^2 u N (P N , s N )u N (P N , S N ) 



(3.7) 



s N 



and 



(3.8) 



The factors (1/N @ ) 1/2 , {l/N T ) 1 / 2 and (l/iVV) 1/2 serve to normalize the wave functions. (See 
the Appendix.) We determine that N N = 0.316, N T = 0.0982, and calculate Nq for each 

— * 

choice of the pentaquark wave function. We write, with k = \k\, 



( 



Vo + 



A 



2{k-k c ) 2 \ 



In 

"W 2 



(3.9) 



and 



2(k-k c ) 2 \ 



Vo + 



\ 



w 



(3.10) 



/ 



We have determined y , A, k c and w from a fit to the wave function given in Fig. 5 of 
Ref.[25]. We find y = -3.66, k c = -0.013 GeV, w = 0.660 GeV and A = 27.73 GeV. For 
the pentaquark we write 

2(k-P N ) 2 

-. 14 

* e (k-P N ) - 



w. 



'Nq I IT 



e 



(3.11) 



where wq is a variable in our analysis. (Note that Nq depends upon the choice of w@.) Here 
k — P N is the relative momentum of the diquark and triquark when P — 0. 
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In the case of a positive-parity pentaquark we assume that the pentaquark decays to a 
positive-parity diquark and a positive-parity triquark. As before, the triquark and diquark 
have zero relative angular momentum. In this case we need to insert factors of 275 at 
the triquark-meson vertex where the quark of momentum k is emitted. In addition, the 
calculation of the normalization factor N T is modified. (See the Appendix.) In Eq. (3.4), 
the trace becomes 



Trace = Tr [(ft - f N + f + m T )ij 5 (lt on + m q )(f N + m N )j° 
(p' on + m q )i l5 (jl' -f N + f + mr) 7 (f + m&)] . 

We may define ^ = 7°^7°, etc. Thus 



(3.12) 



Trace = Tr 



(# - f N + f + m T )^ on + m q ){f N + m N )(^ on + m q ) 



7 5 (^ -P N + P + m T )(P + m e ) 



Tr 



(# - f N + f + m T ){-% n + m q )(-f N + m N )(-k' on + m q ) 
-f N + f + m T )(f + m e ) 



(3.13) 



(3.14) 



In this case we find N N = 0.316 and N T = 0.0201, where only the second value has 
changed relative to the values given above. 



IV. NUMERICAL RESULTS 



In this section we present the results of our numerical calculations. In Fig. 3 we exhibit 
the width for the decay 6 + — > N + K + in the case that the pentaquark has positive parity. 
The result for a negative-parity pentaquark is shown in Fig. 4. The parameters used are 
given in the caption of Fig. 3. The small value for the diquark mass used here corresponds 
to the analysis made in Ref. [29]. In the case of the negative parity pentaquark (Fig. 4) 
very small widths may be obtained for w in the range 0.75 GeV < w < 0.80 GeV. For the 
case of positive parity, quite small widths are found in the range 0.70 GeV < w < 0.90 GeV. 
[See Fig. 3.] 

Results for the decay 0° — > P + D*~ are given in Figs. 5 and 6. Again, we find very 
small widths for the negative-parity pentaquark in the range 0.75 GeV< w < 0.80 GeV with 
widths of about 3^4 MeV for the positive-parity 6° in the corresponding range. (The 
parameters used are given in the figure captions.) 
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FIG. 3: The width for the decay + — > N + K + for a positive-parity pentaquark as calculated in 
our model is shown. Here m e = 1.540 GeV, m T = 0.800 GeV, m D = 0.400 GeV, m N = 0.939 
GeV and m K = 0.495 GeV. 
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FIG. 4: The width of the negative parity pentaquark is shown as a function of the parameter w. 
(The quite small values in the vicinity of the minimum have rather large uncertainties because of 
the limitation of the number of points used in our five-dimensional integral which determines the 
width.) The parameters used are given in the caption to Fig. 3. 
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FIG. 5: The calculated width of the positive-parity charmed pentaquark is shown for the decay 
9° -► P + D*-. Here. m e o=3.099 GeV, m P =0.939 GeV, m D = 0.400 GeV, m T = 2.00 GeV, 
m q = 0.350 GeV and m D *- = 2.007 GeV. 



0.030 -, 




w (GeV) 

FIG. 6: The calculated width of the negative-parity charmed pentaquark is shown. The calculation 
is made for the parameters given in the caption of Fig. 5. 
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FIG. 7: The calculated width for the decay E — > E + it is shown for the case that the S 
has positive parity (dashed line) and for negative parity (solid line). The parameters used are 
mg=1.862 GeV, m T =0.800 MeV, m D = 0.700 GeV, m s = 0.450 GeV, m H = 1.321 GeV and 
m w = 0.139 GeV. 

As stated earlier, the observation of the S pentaquark is a matter of some controversy. 
Its observation was reported in Ref. [3] with a width less than 18 MeV, however, that state 
was not seen in the work reported in Refs. [5] and [11]. Our results calculated for the decay 
S — > 2~ + 7r~ are shown in Fig. 7, while our results for the decay S — > £~ + if - 
are given in Fig. 8. In the case of a negative parity S the summed results for the 
two decay channels are given in Fig. 9 (negative-parity pentaquark) and Fig. 10 (positive- 
parity pentaquark). The calculated widths are of the order of 4-8 MeV in the region 0.75 
GeV< w < 0.80 GeV. The S is not seen in some experiments that involved a search for 
this pentaquark [5, 11]. 

In Fig. 11 we present our results for the decay N + — > N + 7r + . Here we use the mass 
suggested in Ref. [12] for the N + . In this case very small widths are calculated. Much larger 
widths for the N + decay to A + K + are found and the calculated values are exhibited in 
Fig. 12. Since the N + and A s + have the quantum numbers of the nucleon, they are expected 
to mix with neighboring states and be hard to distinguish from other nucleon resonances. 
The mixing of such states is discussed in Ref. [32]. 
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w (GeV) 

FIG. 8: The calculated width for the decay S — ► £~ + K~ is shown for the case that the E~~ 
has positive parity (dashed line) and for negative parity (solid line). The parameters used for mg, 
tut, mo, and m q are given in the caption to Fig. 7. Here m s - = 1.197 GeV and = 0.495 
GeV. 
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FIG. 9: The calculated widths for the decays S — > H~ + 7r~ and S~ + are summed and 
shown in the figure. Here the S is assumed to have positive parity. The parameters used are 
given in the captions of Figs. 7 and 8. 
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FIG. 10: Same caption as Fig. 9 except that the S is taken to have positive parity. 
V. DISCUSSION 

It is of interest to note that the same computer code can yield different values for the 
pentaquark widths depending upon the mass values inserted for the pentaquark, final-state 
baryon and meson, and for the diquark and triquark. We do recognize that our model is 
limited since we have not taken into account the full symmetry of the wave function. That 
is, we do not consider the identity of the quarks in the diquark and the triquark. It is 
possible that a more complex calculation made in the future may overcome this limitation 
of the model. On the other hand, our results are quite suggestive, since very small widths 
are obtained for the 6 + and 6° in accordance with the experimental situation. It is clear 
that further experimental studies and theoretical analysis is desirable. 

APPENDIX A: NORMALIZATION PARAMETER FOR THE NEGATIVE- 
PARITY AND POSITIVE-PARITY PENTAQUARK STATES 

The normalization of the final-state baryon wave function is calculated by requiring that 
that the baryon contain a single quark in addition to the scalar diquark. In that calculation 
the diquark appears on mass shell. We find in the case of a nucleon in the final state, the 
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FIG. 11: The values obtained for the decays iV + — > N + ir + are shown for a positive-parity N + 
(dashed line) and for negative parity (solid line). The parameters used were uim = 0.939 GeV, 
m n = 0.139 GeV, m N + = 1.440 GeV, m T = 0.800 GeV and m d = 0.350 GeV. 




FIG. 12: The calculated width for the decays 7V+ — > A + K + are shown for a positive parity 
(dashed line) and for negative parity (solid line). The parameters used were m^o = 1.116 GeV, 
m Ns = 1.700 GeV, m D = 0.400 GeV, m s = 0.450 GeV, m K = 0.495 GeV and m T = 1.00 GeV. 
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normalization parameter is 

N N = ^-J -^(kk°)(k° + rn q )\* N (k)\ 2 , (Al) 
^ Jo E D {k) 

where E D (k) = (k 2 + m 2 D ) 1 / 2 , with m D being the mass of the diquark and k° = m N — E D (k). 
Here we take k max = 0.7 GeV. In Eq. (Al) 

/ £ _2(k-k c ) 2 \ 

®N{k) =\yo + 7= exp »o . (A2) 



[See Eq. (3.9).] Values of y , w , k c and A are given after Eq. (3.10). 

The calculation of the normalization factor for the pentaquark yields a similar result. For 
example 

N @ = ^-J -^*°*(*° + m T ) | * e (*) | 2 , (A3) 
^ Jo E D (k) 

with k° = rriQ — En(k). Here Eo(k) = (k 2 + m 2 J ) 1 ^ 2 , with mr, being the mass of the diquark. 
For the pentaquark, we use the wave function 

A 2(fc-P JV ) 2 

*e(fc) = 7= exp ~ (A4) 

in Eq. (A3). [See Eq. (3.11).] Here w is a variable of our model. Results for various values 
of w are given in the figures for different pentaquarks. 

The triquark normalization factor is calculated by requiring that the triquark contain 
a single quark in addition to the final-state meson which is placed on mass shell in the 
calculation. We find 

N T = — k 2 dk(k 2 + k 2 + 2k m q + ml) | tf T (fc) | 2 , (A5) 

4n z J 

where k = m T — E K {k). Here E K = (k 2 + m 2 K ) l l 2 , where mx is the mass of the final-state 
meson, and m q is the mass of the exchanged quark. [See Fig. 1.] Note that, in this case, 

j\ _2(k-k c ) 2 

*r(fc) = Vo + 7= exp < , (A6) 

with the parameters yo, wq, k c and A given after Eq. (3.10). 

If we consider positive parity pentaquark states, we do not have to change the baryon 
normalization parameter which was given in Eq. (Al). 
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The normalization parameter of the pentaquark wave function is still of the form given 
in Eq. (A2), while the triquark normalization factor is now 

Nt= 4^J *? d H*v + &~ 2k v m i + m 2 q ) \ 4f T (k) | 2 . (A7) 

Comparison to Eq. (A5) shows that the sign of the term 2fc m g in Eq. (A7) has changed 
relative to the sign in Eq. (A5). 
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